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Abstract 

Given a non-hermitean matrix M, the structure of its minimal poly- 
nomial encodes whether M is diagonalizable or not. This note will ex- 
plain how to determine the minimal polynomial of a matrix without going 
through its characteristic polynomial. The approach is applied to a quan- 
tum mechanical particle moving in a square well under the influence of a 
piece-wise constant PT-symmetric potential. Upon discretizing the con- 
figuration space, the system is decribed by a matrix of dimension three. It 
turns out not to be diagonalizable for a critical strength of the interaction, 
also indicated by the transition of two real into a pair of complex energy 
eigenvalues. The systems develops a three-fold degenerate eigenvalue, and 
two of the three eigenfunctions disappear at this exceptional point, giv- 
ing a difference between the algebraic and geometric multiplicity of the 
eigenvalue equal to two. 

1 Introduction 

Genuinely PT-invariant operators may or may not possess a complete set of 
cigcnstates. In other words, PT-invariance of a matrix M is compatible with the 
presence of (non-trivial) Jordan blocks while hermiticity is not. When consider- 
ing a familiy of PT-invariant operators depending on a parameter, their spectra 
often change qualitatively if one passes through an exceptional point |X] where 
diagonalizability breaks down. It is thus important to be able to either check 
whether a given matrix is diagonalizable, or to locate exceptional points when 
presented with a continuous family of matrices. 

The purpose of this note is to describe a method which allows one to identify 
exceptional points of finite-dimensional non-hermitean matrices by means of an 
algorithm. It is different from the method outlined in [21 EI as it directly aims 
at the minimal polynomial containing the relevant information about (non- 
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) diagonalizability. While being more transparent in the first place, it also 
requires no knowledge of the characteristic polynomial of the given matrix. 

The presentation to follow is problem-based: the algorithm will be developed 
while studying a specific example, the discretized PT-invariant square well. This 
physical system is introduced in Section [3 and it is subjected to the test for 
diagonalizability in the subsequent section. The results will be discussed in 
Section 0] and some open questions will be addressed. 

2 The discretized PT- symmetric square well 

Consider a quantum particle in a one-dimensional box of length 4L subjected 
to a piece-wise constant PT-symmetric potential, 

( -iZ, -2L < x < 0, 
V{x) = l 0, x = 0, ZeR, (1) 

[ iZ, 0<x< 2L, 

which has proved a useful testbed for the discussion of PT-symmetric systems. 
Its eigenvalues are given as the zeros of a transcendental equation jj], and, for 
each value of Z, the lowest two real eigenvalues are known to coalesce and then 
disappear jointly at critical values of the parameter. 

Let us introduce a toy-version of this system by discretizing its configura- 
tion space. This strategy has been applied successfully to decribe tunneling 
phenomena in a driven double-well potential in terms of a three-state model 
5 . Effectively, this technique corresponds to turning Feynman's "derivation" 
of Schr0" dinger's equation from a discrete lattice 6 upside down. Explicitly, 
the continuous set of points of configuration space with labels —2L<x<2L, 
are replaced by five equidistant points at 0, ±L, and ±2L. The wave function 
is allowed to take nonzero values only at these points, so it will be a vector 
with five components at most. However, the hard walls of the square well at 
±L force the wave function to vanish there, leaving us with only three non-zero 
components, tp — ► tpk = ip(kL),k — 0, ±1. The potential energy defined in (f3~|) 
turns into a diagonal matrix, V(x) — ► V =diag(— iZ, 0, iZ). The operator for the 
kinetic energy follows from replacing d 2 ip(x)/dx 2 — > (ipk+i — %4>k + tpk-i)/L 2 . 
Putting all this together, the Hamiltonian operator of the discrete version of 
this system reads 

/ 1 \ 

H ~ 2E - H , where H = 1 1 , f = Z/t] ; (2) 

V 1 -%£ ) 

here E is the (3 x 3) identity matrix, and an overall factor -q — h 2 /2mL 2 has 
been dropped. The matrix H inherits PT-invariance from the square well: the 
matrix H, and hence Ho, is invariant under the combined action of parity P, 
represented by a matrix with unit entries equal to one along its minor diagonal 
and zero elsewhere, and T, effecting complex conjugation. In the next Section, 
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the diagonalizability of H, the nontrivial part of the Hamiltonian Ho, will be 
studied. 



3 Diagonalizability of the PT-symmetric square 
well 

The minimal polynomial tom of a matrix M is defined (see for example) as 
the polynomial of least degree in M which annihilates M, that is, tom(M) = 0. 
This polynomial is unique if the coefficient of its highest power is taken to be 
one: the minimal polynomial is monic. Since any matrix M of size N, say, is 
annihilated by its own characteristic polynomial, pm(M) = ; the degree of the 
minimal polynomial does not exceed N. Once Tom has been found, one needs 
to determine whether it has only single roots, i.e. whether 



holds. If it does, the matrix M is diagonalizable - otherwise, it is not diagonal- 
izable since multiple roots of mM indicate the presence of Jordan blocks larger 
than one. 

The procedure to determine the minimal polynomial of M, as outlined in 
invokes the characteristic polynomial of M and repeated applications of the 
Euclidean algorithm generalized to polynomials. The method presented below, 
taken from jS] , aims at directly constructing the minimal polynomial. The fun- 
damental observation is that matrices of size N constitute a vector space of 
dimension iV 2 . This is seen immediatly by setting up a one-to-one correspon- 
dence between matrices of size 3 and vectors of length 9 = 3 2 , for example, 
simply by rearranging the elements Mjk of each M systematically according to 



In view of this correspondence, the theorem by Cayley-Hamilton - every matrix 
satisfies its own characteristic equation, pm(M) = 0, - turns into a statement 
about linear dependence of the (N + 1) vectors E = M°, M, M 2 , . . . , M*. In 
order to find the minimal polynomial of a matrix M one thus simply calculates 
M™, n = 1 . . . N, and then determines successively whether the vectors E and M 
are linearly independent; if not, one adds M 2 and asks the same question; etc. 
Proceeding in this way, one is obviously able to identify linear dependence among 
the vectors M™ containing only the smallest powers necessary. This, however, 
comes down to the definition of the minimal polynomial of M. Gram-Schmidt 
orthonormalization effectively provides a systematic test for linear dependence 
among the first k elements of M n ,n = . . .N. 



uo(<N) 

rn M (A) = J| (A — M v ) , all M v distinct . 



(3) 



M <^> (Mn, M12, M13; M 2 i, M 22 , M 23 ; M 3 i, M 32 , M 33 ) T . 



(4) 
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Applying these ideas explicitly to the matrix H in J3J) leads to 

E & (1,0,0;0,1;0; 0,0,1), (5) 

H & (if, 1,0; 1,0,1; 0,1, -if), (6) 

H 2 «• (l-f 2 ,-if,l;-i£,2,if;l,if,l-f 2 ), (7) 

H 3 «*■ (2-f 2 )(if,l,0;l,0,l;0,l,-if) =(2-f 2 )H. (8) 

It is easy to see that neither the first two nor the first three vectors in this 
sequence are linearly dependent. Consequently, there must be a relation ex- 
pressing H 3 in terms of the others, and indeed, the minimal polynomial follows 
immediately from 

(2-£ 2 )H-H 3 = =>. mH (A) = A 3 + (£ 2 - 2)A . (9) 

In addition, the characteristic polynomial of H must coincide with wh since it 
there is only one monic polynomial of third degree annihilating H. 

If one is not able to factor the resulting minimal polynomial, one needs 
to check whether the minimal polynomial and its derivative m'i-i(A) have a 
common factor which can be achieved by applying the Euclidean algorithm to 
this pair (cf. In this present case, this amounts to writing mn(A) = 

(A - a)m' H (\)/3 + Ri(X), implying that a = and Ri(X) = (2/3) (f 2 - 2)(A + 
1/2). Two different cases arise: if f ^ 2, one finds that the only common 
factor of the minimal polynomial and the derivative is equal to one - thus, the 
minimal polynomial is of the form and the matrix H must have three different 
eigenvalues making it diagonalizable. If f 2 = 2, the algorithm immediately stops 
and thus identifies A 2 as the highest common factor of mn(A) and its derivative. 
Consequently, the minimal polynomial has a three-fold root A(= 0), indicating 
that H is not diagonalizable. 



4 Discussion and Outlook 

The properties of H at the exceptional points defined by f± = ±v2 deserve 
a brief discussion. It is not difficult to see that the geometric multiplicity of 
the eigenvalue is one at the exceptional points (H has only one non-zero 
eigenvector) while its algebraic multiplicity equals three (zero is a triple root 
of toh). Contrary to previously studied cases, three eigenvalues coalesce for 
£ — * f±. For all values of f , the vector (1, —if, —1) is an eigenstate of H(f) with 
eigenvalue 0, and it does not exhibit any particular behaviour at the critical 
values f±. Therefore, it seems reasonable to say that the eigenstates associated 
with the two A-dependent eigenvalues disappear at the exceptional point. It is 
not obvious from a physical point of view why this scenario is preferred over the 
familiar situation of just one disappearing eigenstate. 

The natural question to ask now is whether one can expect algorithmic 
tests for diagonalizability to exist for a quantum system living in a Hilbert 
space accomodating a countable infinity of states. As one needs to potentially 
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perform an infinite number of steps, the idea of a useful algorithm gets somewhat 
blurred. Nevertheless, it is likely that one can search for systematic properties 
of finite-dimensional approximations which, hopefully, behave smoothly in the 
limit of infinite dimension. 
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